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ABSTRACT 



We are conducting a high angular resolution imaging survey of Galactic Wolf- 
Rayet stars using the Wide Field and Planetary Camera 2 aboard the Hubble 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. 

2 Present address: NASA/GSFC, Code 685, Building 21, Greenbelt, MD, 20771, wal- 
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Space Telescope. We have found a small stellar cluster associated with the faint, 
close pair WR 38 and WR 38a. We present astrometric measurements and pho- 
tometry in the wide-band F336W (U), F439W (B), and F555W (V) filter system 
for these cluster and nearby stars. We combine our photometry with Johnson 
and IR magnitudes, and compare the observations with calibrated model results 
for reddened stars to adjust the HST zeropoints and to identify five probable 
main sequence members of the cluster. A least-squares fit of the colors and mag- 
nitudes of this set yields a cluster reddening of E(B — V) = 1.63 ± 0.05 mag 
and a distance of lO^ 2 kpc for an assumed ratio of total-to-selective extinction 
of R = 3.1. We discuss the relationship of this cluster to other objects along the 
line of sight. If situated at a distance of ~ 8 kpc, then the cluster would reside 
in a dense region of the Carina spiral arm, close to a giant molecular cloud and 
the starburst cluster NGC 3603. 

Subject headings: binaries: visual — stars: imaging — stars: Wolf-Rayet - 
stars: individual (WR 38, WR 38a) — open clusters and associations: individual 
(C1104-610a) 

1. Introduction 

Wolf-Rayet (WR) stars are the evolved, He-burning cores and hot stellar envelopes of 
once-massive O-type stars. Their strong winds (terminal velocities: Voo ~ 1000 - 2500 km s _1 ) 
and heavy mass loss rates (M — ~ [1-6] x 10~ 5 M yr" 1 ; Willis 1991; Nugis & Lamers 2002) 
power their characteristic strong, broad emission lines of He, N, C, and O, which correspond 
to increasingly evolved evolutionary states. While the environment and multiplicity of O- 
type stars have been explored in depth (Gies 1987; Mason et al. 1998), these characteristics 
for the evolved WR phase are not as well determined. In 1996 we initiated a survey with 
the Hubble Space Telescope Wide Field and Planetary Camera 2 (HST WFPC2) to discover 
and quantify the multiplicity and environments of Galactic WR stars. Here we discuss two 
of our targets, WR 38 and WR 38a, that were revealed to be both multiple and contained 
within a small stellar cluster (Wallace et al. 1999). 

WR 38 is a WC4-type star (MacConnell & Sanduleak 1970) and WR 38a is a WN5-type 
(Shara et al. 1991). The two are separated by only 20" in the sky, and they are neighbors to 
another even tighter 5" pair of WR stars, WR 38b and WR 39, at a separation of 3'8. These 
four plus WR 37 are surrounded by a ring nebula (Marston 1997), although it is not clearly 
associated with any one of these stars. Our view in this direction in the Galactic plane 
cuts through the Sagittarius-Carina spiral arm at distances of approximately 2.4 and 8 kpc 
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(Grabelsky et al. 1988; Georgelin et al. 2000; Russeil 2003). The photometric distances of 
these WR stars are generally close to 5 kpc (van der Hucht 2001), which would place them 
among the young population of the Carina spiral arm. However, in a recent paper Shorlin, 
Turner, & Pedreros (2004) presented photometry for the stars in the immediate vicinity of 
WR 38 and WR 38a that indicated a much greater distance of 14.5 ±1.6 kpc. This distance 
would place the pair in an extension of the outer Perseus arm some 14 kpc from the Galactic 
center. There are only 4 of 227 known Galactic WR stars that have such a great distance 
from the Galactic center (van der Hucht 2001). Shorlin et al. (2004) designate the cluster as 
C1104-610a, and images of the field surrounding the cluster are presented in Wramdemark 
(1976) (Plate 4), van der Hucht et al. (1981), and Shorlin et al. (2004) (Fig. 1 and 2). 

Here we present new astrometric and photometric results on the stars found close to 
WR 38 and WR 38a based upon HST 'PC images with superb angular resolution. We identify 
which stars are associated with the WR pair based upon their derived reddening and location 
in a color-magnitude diagram. Our results are mainly consistent with those of Shorlin et al. 
(2004) but we argue that the acceptable distance range is such that the WR pair and its 
surrounding cluster may reside in the distant portion of the Carina spiral arm. 

2. Observations and Data Reduction 

Observations of the fields of WR 38 and WR 38a were obtained in 1996 March (JD 
2450166.58) through the F336W, F439W, and F555W filters as part of our HST WFPC2 
survey of Galactic WR stars. These WFPC2 filters were designed to correspond to the 
Johnson UBV filter set (Johnson & Morgan 1953; Johnson et al. 1966), but because of the 
strong emission lines present in WR spectra the transformation from the HST to the Johnson 
system is approximate at best for them. Additionally, the WFPC2 F336W filter suffers from 
a red leak which may cause dusty WR or other very red stars to appear brighter than in a 
similar U band image. Therefore, all of the data reductions and analyses were performed 
within the WFPC2 synthetic magnitude system (Holtzman et al. 1995). 

We centered the WR stars within the PC chip to take advantage of the high resolution 
(scale of 0'.'046 pixel -1 ). The observations also utilized a high gain setting, with a threshold 
of ~ 53000 photoelectrons pixel -1 , to optimize the discovery of faint companions. For each 
observation we made two exposures to facilitate the removal of cosmic rays (40 s for the 
F336W, 10 s for the F439W, and 2 s for the F555W filter). Since these are SNAPSHOT- 
type observations, which are designed to fill holes in the HS T schedule, the observations were 
taken using only one guide star. Fortunately, due to the short exposure times, this resulted 
in no discernible loss of data quality due to telescope motion, although this did place limits 
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on the accuracy of the absolute astrometry as the telescope may roll during the observation 
(Baggett 2002; Gonzaga 2002; Nelan & Makidon 2002). 

WR 38 is sufficiently close to WR 38a that both stars could be captured within the same 
WFPC2 image (Fig. 1). The HST image clearly shows a number of components surrounding 
both WR stars that are not fully resolved in ground-based images (Shara et al. 1991; Shorlin 
et al. 2004). Without spectroscopic data we can only estimate the identity of the WR 
star within each stellar grouping by assuming it is the brightest object within the existing 
astrometric error box. The numbering system we adopted is based upon angular separation 
from WR 38 (star #1), and WR 38a is labeled as #21. 

After normal pipeline processing (Baggett 2002), we aligned the six raw images using 
the task IM ALIGN in IRAF 1 /STSDAS. We then combined each pair of images per filter 
using the routine CRREJ. We identified the stars in each combined image using the task 
DAOFIND with the detection threshold set at 10<r over the background level. We determined 
a background value by using the IMEXAMINE routine to derive the mean background value 
within 5x5 pixel boxes located at 5-10 well-distributed positions across the image. The 
single average of these values was used for the detection threshold. 

Because the observations were made with just one guide star, the images may suffer 
from telescope roll and jitter. Consequently, the point spread function (PSF) of the stellar 
images may be asymmetrical, leading the DAOFIND algorithm to select false stars. Thus, 
all the detected stars were visually inspected before acceptance. 

We used this list as input to the IRAF/STSDAS package METRIC to derive the as- 
trometric information given in Table 1. METRIC translates the WFPC2 pixel positions 
to celestial coordinates after correcting for geometric distortion. The relative positions of 
well-exposed stars are accurate to better than 0'.'005 for targets imaged on one chip (er- 
rors are primarily due to centering uncertainties). However, while the relative positions 
are highly accurate, the absolute positions of the objects are subject to larger errors of 
±(0.5 — 5) arcseconds due to the use of a single guide star (Baggett 2002; Gonzaga 2002; 
Nelan & Makidon 2002). Fifteen of the stars in Table 1 appear in the 2MASS All-Sky Cat- 
alog of Point Sources (Cutri et al. 2003), and the mean differences in the coordinates are 
a(2MASS) - a(HST) = +0.175 ± 0.008 s and 5(2MASS) - 5(HST) = +0"14 ± 0"18. 

We performed aperture photometry using the IRAF package DAOPHOT (Stetson 1987). 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 
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We chose aperture photometry over PSF photometry because of the difficulties involved in 
fitting the PSF due to the small telescope motions. We used an 11 pixel (0'.'5) radius aperture 
to determine the total stellar flux, and we estimated the sky background using the OFILTER 
algorithm within the PHOT package of DAOPHOT (as recommended by Ferguson 1996). 

A number of photometric corrections are required to transform from raw to calibrated 
magnitudes (Baggett 2002). The stellar magnitudes were corrected for charge transfer ef- 
ficiency using the Whitmore, Heyer, & Casertano (1999) formula for faint sources and the 
Dolphin (2000) 2 formula for more luminous sources (brighter than 100e~). Geometric dis- 
tortion effects were corrected using a scheme that rescaled the flux measurements of the 
stars based on their positions within the image. Contamination corrections were made using 
information taken from the WFPC2 Web site 3 . The 34 th row effect was corrected using the 
scheme of Anderson & King (1999). Most significantly, the zero-point corrections for the 
F555W and F439W magnitude data were calculated using the Dolphin (2000) updates to the 
original Holtzman et al. (1995) offsets. The Dolphin (2000) revisions to the Holtzman et al. 
(1995) zero-points (~ 0.01 mag) are not the result of temporal effects, but rather of improved 
calibration methods, and hence, they have greater accuracy. Unfortunately, Dolphin did not 
provide updates for the F336W filter, so we used the Holtzman et al. (1995) calibration for 
the F336W filter images. 

Our final magnitudes in the HST system are collected in Table 2. We have included 
a transformation of the HST wide-band magnitudes to Johnson BV using the scheme of 
Dolphin (2000). We caution that this scheme may not be reliable for the WR stars since 
both the B and V band-passes are sensitive to the strong emission lines present in WR 
spectra (van Genderen, Steemers, & van der Hucht 1987). We did not attempt to transform 
our F336W magnitudes to Johnson U since the latter straddles the Balmer jump while the 
former is sensitive mainly to shorter wavelength flux. We have also included in Table 2 the 
corresponding identification number (STP) and magnitude information from the work of 
Shorlin et al. (2004) for 7 stars in common (mainly cluster members). 

We found that there are systematic differences between our transformed V and B — V 
magnitudes compared to those of Shorlin et al. (2004). We omitted the WR stars from 
this comparison because of our concerns about the transformation to Johnson magnitudes, 
and this leaves 5 stars in common between our samples. We find mean differences of < 
Vust - V ST p >= 0.24 ± 0.03 and < (B - V) U st - (B - V) STP >= -0.12 ± 0.09, i.e., the 
mean differences are larger than the standard deviations from the mean differences (and the 



2 http: / /www. noao.edu/stafF/dolphin/wfpc2_calib/ 
3 http: / / www.stsci.edu/instruments/wfpc2 / 
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standard deviations are comparable to the expected observational errors). These differences 
are probably due to errors in the adopted zero-points for the WFPC2 images. We have 
found, for example, that there is an intrinsic scatter of order 0.1 mag between the observed 
and synthetic HST magnitudes for some two dozen WR stars with NUV spectroscopy from 
the International Ultraviolet Explorer Satellite and optical spectrophotometry (Torres & 
Massey 1987; Torres-Dodgen & Massey 1988). While part of this scatter may be due to 
errors in the calibration of the spectrophotometry, we caution that zero-point offsets of order 
0.1 mag may be required to place our photometric results on an absolute scale. Thus, the 
relative magnitudes for stars within a given image are secure, but color indices formed from 
magnitude differences between images may require a zero-point shift. In the next section we 
describe one way of estimating the zero-point adjustments based upon a reddening analysis 
of multi-wavelength photometry. 



3. Interstellar Reddening 

Extinction due to interstellar dust decreases at longer wavelengths, so the best approach 
to determining reddening and extinction is to compare both long and short wavelength flux 
measurements with spectral models transformed using a suitable reddening law. We were 
pleased to learn that 15 of the 25 targets in our HST field are included in the 2MASS All-Sky 
Catalog of Point Sources (Cutri et al. 2003) that lists J, H, and K magnitudes for most of 
them. Here we examine the 2MASS, HST, and, where available, Johnson magnitudes from 
Shorlin et al. (2004) in order to estimate the reddening for each target. Our approach is to 
calculate synthetic magnitudes in each band based on a consistent spectral flux distribution 
for a given ZAMS stellar effective temperature, T eff , and a given reddening, E(B — V). 

We assembled a small grid of theoretical spectra covering the wavelength range of interest 
for ZAMS stars along the upper main sequence. We used model spectra based on non-LTE, 
line-blanketed model atmospheres from Lanz & Hubeny (2003) for the O-stars (for effective 
temperatures of 27500, 35000, and 45000 K), and we adopted spectra from LTE models 
by R. L. Kurucz 4 for the B-stars (for 10000, 15000, 20000, and 24000 K). The gravity 
log g for each temperature was taken from the 1 Myr ZAMS sequence of Lejeune & Schaerer 
(2001). In addition to the computational differences between the Kurucz and Lanz & Hubeny 
models, the former assume a microturbulence of 4 km s^ 1 while the latter use 10 km s _1 . 
Despite these differences, the predicted flux distributions generally make a smooth transition 
at the temperature boundary between them (where non-LTE effects for hydrogen become 



4 http://kurucz. harvard.edu/grids/gridP00/fp00k4.pck 
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negligible). 

We calculated synthetic magnitudes for these spectra using the SYNPHOT routine 
CALCPHOT (Bushouse & Simon 1998), which is distributed with the STSDAS package 
in IRAF. The software includes the spectral response functions for the HST wide band 
filters as well as those for the standard Johnson filters. However, we used the response 
functions from Cohen, Wheaton, & Megeath (2003) to simulate the sampling of the spectra 
in the IR with the 2MASS filter system. The procedure returns synthetic magnitudes on 
the VEGAMAG system, so defined that Vega has a magnitude 0.0 in all filters. Since 
Vega actually has magnitudes U = 0.02, B = 0.02, and V = 0.03 in the Johnson system, 
we adjusted the both the Johnson and corresponding HST synthetic magnitudes by these 
amounts to make them consistent with observations. The resulting synthetic ZAMS colors 
relative to Johnson V and the absolute K magnitude are collected in Table 3 for our grid 
of stellar effective temperatures. These are generally in reasonable agreement with Johnson 
color indices published by Wegner (1994) and Bessell, Castelli, & Plez (1998) except perhaps 
for the shortest wavelengths and hottest stars. For example, Turner (1976) and Underhill 
et al. (1979) argue that the hottest O-stars have B — V — —0.32 versus our lower limit of 
B — V = —0.28. Note that the U — V and B — V color indices are often based on uncertain 
de-reddening procedures for O-type stars, which may tend to produce colors that are too 
blue. The lowest B — V = —0.28 for the 45000 K model is consistent with that for the bluest 
O-stars in the Galaxy (Maiz-Apellaniz et al. 2004) and in the LMC (Fitzpatrick 1988). 

We next transformed each of these model spectra for a grid of 40 values of E(B — 
V) from to 3.9 using the reddening scheme of Fitzpatrick (1999) for a ratio of total-to- 
selective extinction of R = 3.1 (the same as adopted by Shorlin et al. 2004). We then used 
CALCPHOT to determine the predicted color indices as a function of the reddening. A set 
of sample reddening curves is shown in Figure 2 for the 35000 K model spectrum. Note that 
unlike the other filters the color index trend for the HST F336W filter undergoes a reversal 
at large reddening because of a red-leak in this filter. 

We performed a numerical grid search to find the stellar temperature and reddening 
that best matched a given set of IR and optical color indices. We found this scheme was 
able to confirm the known reddening in a number of test cases for OB stars where Johnson 
U, B, and V and 2MASS J, H, and K magnitudes were available, but the method is not 
applicable to stars with temperatures cooler than our grid limit of 10000 K. We applied the 
method to all the stars in our WFPC2 field for which 2M ASS photometric data exist. These 
stars are listed in Table 4, which gives our assigned number, the 2MASS designation, the 
number adopted by Shorlin et al. (2004), our derived reddening, and de- reddened estimates 
of m(336) — K and K (see below). There are several stars with relatively low reddening 
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that are probably foreground objects, but most of the other stars have a reddening in the 
range E(B — V) = 1.5 to 2.0 (with the possible exception of the faint star #8). There are 
four cluster stars with a complete set of WFPC2, Johnson (Shorlin et al. 2004), and 2MASS 
magnitudes (after omission of the two WR stars for which the method may be poor), and 
the mean reddening for these four stars is E(B — V) — 1.63 ± 0.05 (mean error), in good 
agreement with the result of E(B — V) — 1.60 ± 0.02 from Shorlin et al. (2004). 

We can use these reddening fits to return to the question of magnitude zero-point off- 
sets that we raised in the previous section. Table 5 lists the average values of the residuals 
from the fit for each filter based on the same sample of four stars with complete magnitude 
coverage. These represent any systematical offsets between the observed magnitudes and the 
synthetic magnitudes derived from the theoretical spectra, so they encompass any lingering 
problems in the observational zero-points, filter responses, and wavelength specific deficien- 
cies in the theoretical models. The offsets in Table 5 are generally small and comparable to 
the scatter in the residuals, but there are some interesting exceptions. We see, for example, 
that the mean difference in the residuals for m(555) — V is 0.21 ± 0.04, which is the same 
within errors as the difference between the transformed Vhst and Vstp that we found in the 
last section. Similarly the residuals in m(439) — m(555) suggest that our uncorrected color 
index is too blue in the same way as demonstrated in our comparison of our transformed 
B — V with the colors from Shorlin et al. (2004). Thus, the mean residuals in Table 5 lead 
us to revise our original WFPC2 zeropoints to 18.305 for the F336W filter, 20.036 for the 
F439W filter, and 21.496 for the F555W filter. An application of these zero-point corrections 
should provide consistent results across for the entire set of photometric measurements. We 
give in columns 5 and 6 of Table 4 the individually de-reddened color index m(336) — K 
(including a zero-point offset of —0.35 ± 0.09 from the mean residuals in Table 5) and ex- 
tinction corrected K (using a zero-point offset of +0.15 ± 0.07 from Table 5) that we will 
use in the next section to estimate the cluster distance. 

We caution that these results are sensitive to the adopted value of total-to-selective 
extinction R. This parameter is often estimated by comparing the reddening EiB — V) to 
the extinction in the IR, and we used equations A3, A4, and A5 from Fitzpatrick (1999) 
to estimate R for the four stars with complete Johnson and 2MASS photometry. The 
mean derived value, R = 3.25 ± 0.15, is consistent within errors with our assumed value 
of R = 3.1, but the application of reddening curves for R = 3.25 leads to a mean reddening 
of E(B — V) — 1.51 ±0.06 for the same sample of four stars using the scheme outlined above. 
We discuss below how such a revision in the reddening law influences the distance estimate. 
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4. Distance to the Cluster 

We can estimate the distance to the cluster surrounding WR 38 and WR 38a based 
upon the positions of the stars in a de-reddened color and extinction corrected magnitude 
diagram. The spectral flux distributions of massive stars appear similar to the Rayleigh- 
Jeans tail of a blackbody spectrum in the optical and near-IR spectral range, and we need 
spectral line diagnostics (primarily through spectral classification) in order to estimate the 
stellar effective temperature reliably (and thus establish the luminosity of the star relative 
to the main sequence). In the absence of spectral data, the best approach is to use a color 
index based upon the widest wavelength range possible as a temperature parameter. We 
selected the m(336) — K index for this purpose, which spans a range of nearly 3 magnitudes 
over the OB star temperature range (compared with a range of only 0.26 mag in B — V; 
see Table 3). We also decided to use K for the magnitude ordinate in the color-magnitude 
diagram since it is least affected by uncertainties in extinction. 

The color-magnitude diagram for the WR 38 cluster in this magnitude system is illus- 
trated in Figure 3. The m(336)— K colors were de-reddened and the K magnitudes extinction 
corrected using the individual E(B — V) values determined for each star in the last section 
(Table 4). In addition, the zero-point corrections found from the mean residuals to the red- 
dening solutions (Table 5) were applied to bring the observed magnitudes into the synthetic 
magnitude system needed to compare the observed and theoretical color-magnitude dia- 
grams. There are five stars that have similar positions in Figure 3 that probably correspond 
to the cluster main sequence. One bright star, #7, is located above the main sequence, and 
it is either an evolved, luminous star within the cluster or a foreground object. The two WR 
stars are found well to the right of the main sequence, which is probably due to the strong 
emission lines found in the X-band (Figer, McLean, & Najarro 1997). 

We adopted the 1 Myr ZAMS from the models of Lejeune & Schaerer (2001) to determine 
the distance modulus from the difference of the observed and theoretical absolute magnitudes. 
We compared our synthetic color index m(336) — K values from SYNPHOT/CALCPHOT 
for our temperature grid with the corresponding values from the tables of Lejeune & Schaerer 
(2001), and we found a mean difference of 0.04 ±0.06 mag. We added this small offset to the 
model ZAMS color index from Lejeune & Schaerer (2001) in order to intercompare the colors 
consistently in the synthetic system. We then compared the observed and model colors for 
a grid of assumed distance modulus values to determine the best fit estimate. The derived 
distance modulus is 5 \ogd — 5 = 15.01}'q based on the x 2 residuals of the fit (corresponding 
to a distance of d = lOl^ 2 kpc). At this distance, the five main sequence stars in Figure 3 
have colors and magnitudes of ZAMS stars in the temperature range 33000 to 38000 K, 
which are characteristic of O-type stars. 
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Our distance modulus estimate agrees within errors with that derived by Shorlin et al. 
(2004), 15.80 ± 0.25. The small difference in distance modulus between these investigations 
is probably due to differences in the adopted ZAMS isochrone. The ZAMS relation adopted 
by Shorlin et al. (2004) (from Turner 1976) is an empirical relation based upon a somewhat 
evolved population of stars, and it is approximately 0.8 mag brighter at B — V = —0.30 
than the theoretical 1 Myr ZAMS that we adopted from Lejeune & Schaerer (2001). Note 
that Schonberner & Harmanec (1995) verified that the Geneva models used by Lejeune & 
Schaerer (2001) yield absolute magnitudes that agree with those determined for eclipsing 
binaries after a correction is made for evolution from the ZAMS. We also show in Figure 3 
the isochrone for an age of 8 Myr, approximately equal to the maximum age of clusters 
containing WR stars (Massey, DeGioia-Eastwood, & Waterhouse 2001), and a comparison 
of the isochrones shows that we may be underestimating the distance modulus by ~ 1 mag 
if the cluster is signftpificantly older than assumed. Unfortunately, we do not know the age 
of the cluster, and consequently the amount of evolutionary brightening at the upper end of 
the main sequence is unknown. Deeper photometric observations could resolve this problem 
through ZAMS fitting of less massive, unevolved stars. 

If the ratio of total-to-selective extinction is revised upwards to R = 3.25 (§3), then the 
derived reddening is smaller and the intrinsic m(336) — K colors are larger than shown in 
Figure 3. Such a revision would lead to a decrease in distance modulus to a value of 13.5 mag, 
slightly below the 1 a error limit quoted above for the R = 3.1 model fit. Thus, given the 
uncertainties in the cluster age and the reddening law, the range in acceptable distance 
modulus value may be somewhat larger than given above. Taken at face value, our analysis 
does not improve upon the the distance estimate from Shorlin et al. (2004), 14.5 ±1.6 kpc, 
but we caution that their error refers only to internal sources and at this stage it is prudent to 
acknowledge an error range that accounts for systematical errors resulting from assumptions 
about the cluster isochrone and reddening law. 

Our line of sight towards WR 38 (/ = 290?57, b = — 0?92) passes twice through the 
Carina spiral arm (see Fig. 4 of Grabelsky et al. 1988) and we can use information from 
CO radio emission maps to help place the giant molecular clouds and associated massive 
star forming regions along this line of sight. Grabelsky et al. (1988) and Dame, Hartmann, 
& Thaddeus (2001) show that there is a large and rare hole in the near side of the Carina 
spiral arm in this direction that might allow us to see some very distant clusters. This first 
cut through the arm occurs at distance of « 2.4 kpc, and Shorlin et al. (2004) discuss some 
of the sparsely populated field stars they find at this distance. The next molecular cloud 
down the line is found at a distance of 6.8 kpc in the direction I = 290?2, b = — 0?2 (with a 
radial velocity of —1 km s _1 and designated #14 in Fig. 2 of Grabelsky et al. 1988), which 
is located about 1?0 away from the WR 38 cluster. This distance is consistent with the 
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distance error range for the WR 38 cluster and the other nearby WR stars. However, just 
beyond this at 7.9 kpc we encounter one of the largest molecular clouds in the Carina spiral 
arm (with a radial velocity of +22 km s _1 and given as #13 in Fig. 2 of Grabelsky et al. 
1988) where our line of sight crosses the spiral arm a second time. The WR 38 cluster lies 
only 0?4 away from the edge of this cloud. There is very little CO emission in the (I, v) 
diagram corresponding to distances beyond this, so it is doubtful that there are favorable 
environments for massive star formation beyond a distance of 8 kpc in this direction (see also 
Fig. 5 in Russeil 2003). However, McClure-Grimths et al. (2004) present deep radio maps of 
21 cm emission from neutral hydrogen in this region that suggest that some weak emission is 
present from an extension of the Perseus arm (at a distance of ~ 15 kpc and a radial velocity 
of +50 km s -1 ) and a possible outer arm (at a distance of ~ 21 kpc and a radial velocity of 
+ 120 km s _1 ). The outer disks of large spiral galaxies occasionally show clear H I spiral arms 
extending to radii in excess of 30 kpc, while the stellar arms stop at radii of about 10 kpc. 
A classic example is M83, where the H I spiral arms extend three times as far as the stellar 
arms (Tilanus & Allen 1993). Note that there are very few WR stars known in the Milky 
Way Galaxy beyond the solar circle (see Fig. 6 in van der Hucht 2001), where they would 
be easy to discover if they existed there (due to generally low reddening). Furthermore, the 
number ratio of WR to O-type stars decreases dramatically in low metallicity environments 
such as the outer Galaxy, so even if there are fair numbers of O stars in the outer Galaxy, 
there will be relatively many fewer WR stars there. 

Georgelin et al. (2000) discuss the stars that are associated with H II regions in this 
general direction, and they find that the extinction ranges from Ay ps 1.8 at 2.8 kpc to 
Ay ~ 2.6 at 4.2 kpc, and up to Ay ps 4 at 8 kpc. Our derived extinction of Ay = 5.0 ± 0.2 
would tend to suggest that the WR 38 cluster also resides at a distance of ~ 8 kpc or larger. 
The 8 kpc distance would place the WR 38 cluster approximately midway between molecular 
cloud #13 (and associated H II complex 289.3 — 0.6; Georgelin et al. 2000) and molecular 
cloud #17 (H II complex 291.6 — 0.7) which is associated with the starburst cluster NGC 3603 
(Moffat et al. 2002). Note that Sung & Bessell (2004) find a reddening of E(B -V) = l.% 
in the outskirts of NGC 3603 that is comparable to that for the WR 38 cluster. Georgelin 
et al. (2000) identify an H II region 290.487 - 0.814 within the larger 289.3 - 0.6 complex 
that has as exciting stars #48 and #52 from the list of Wramdemark (1976). The two stars 
have a comparable extinction {Ay ps 3.7) and are located only ps 0?3 away from the cluster, 
so these stars and the H II region could be related to the WR 38/WR 38a cluster. 

The solution to the distance problem for the WR 38/WR 38a cluster will come when 
spectra can be obtained for the cluster stars. Classification dispersion spectra would help 
establish the temperatures and luminosities of the stars, and thus determine accurately the 
distance modulus of the cluster. Furthermore, radial velocities from spectra would also allow 
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us to distinguish between the differential Galactic rotation expected for the locations along 
the Carina and Perseus spiral arms. Such observations would clearly help us understand the 
history of massive star formation in this distant region of the Galaxy. 
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Fig. I.— The HST/PC image of WR 38 (#1) and WR 38a (#21) made with the F555W 
filter (placed according to J2000 coordinates). The identification number above each star 
corresponds to the system used in Tables 1 and 2. 



Fig. 2. — Reddening curves based on synthetic magnitudes from SYNPHOT/CALCPHOT 
for a spectrum from a model with T c g = 35000 K and logg = 4.16 (Lanz & Hubeny 2003) 
and for a reddening law with R = 3.1 (Fitzpatrick 1999). The top panel shows how the HST 
color indices m(336) — V, m(439) — V, and m(555) — V vary with reddening E(B — V). The 
m(336) — V trend reverses at large reddening due to the red leak in the the F336W filter. 
The middle panel shows the variations for the Johnson U — V and B — V indices, while the 
lower panel shows the same for the IR indices, J — V, H — V, and K — V. 



Fig. 3. — A color- magnitude diagram based upon the de-reddened color index (m(336) — K) 
and the absolute K magnitude (for a distance modulus of 51ogrf — 5 = 15.0). The solid 
(dashed) line represents the 1 Myr (8 Myr) isochrone from Lejeune & Schaerer (2001) while 
the plus signs represent the magnitudes and colors of the probable cluster stars. Errors in 
the color index are represented by the horizontal line through each symbol (errors in K are 
comparable to the vertical symbol size). The two WR stars that appear in the right hand 
side of the diagram and the bright star (#7) located at the top of the diagram were omitted 
in the main sequence fitting procedure. 



Table 1. Astrometric Data for the WR 38/ WR 38a Cluster 11 



Field 


a 




5 


Sep. b 


P.A. b 


Number 


(J2000) 


(J2000) 




(°) 


1 


11:05:46.62 


-61: 


:13:48.8 






2 


11:05:46.49 


-61: 


43:50.2 


1.63 


215.2 


3 


11:05:46.30 


-61: 


43:45.9 


3.75 


322.1 


4 


11:05:46.74 


-61: 


43:43.3 


5.57 


8.9 


5 


11:05:47.56 


-61: 


13:55.3 


9.38 


133.4 


6 


11:05:47.32 


-61: 


43:40.3 


9.90 


30.7 


7 


11:05:47.36 


-61: 


43:40.1 


10.28 


31.5 


8 


11:05:48.23 


-61: 


43:47.1 


11.79 


81.7 


9 


11:05:46.81 


-61: 


44:00.0 


11.87 


173.2 


10 


11:05:48.24 


-61: 


43:52.8 


12.34 


108.8 


11 


11:05:46.68 


-61: 


43:35.1 


13.69 


1.8 


12 


11:05:45.79 


-61: 


43:35.8 


14.32 


335.3 


13 


11:05:45.33 


— Ol. 


•1 ^-^fi Q 
lo.ou.y 


15.08 


322.0 


14 


11:05:48.63 


-61: 


43:41.7 


16.21 


63.8 


15 


11:05:48.87 


-61: 


43:46.8 


16.34 


82.8 


16 


11:05:46.41 


-61: 


13:31.3 


17.59 


355.1 


17 


11:05:46.21 


-61: 


43:31.1 


17.95 


350.6 


18 


11:05:47.72 


-61: 


13:32.5 


18.17 


26.0 


19 


11:05:46.78 


-61: 


43:30.3 


18.58 


3.6 


20 


11:05:44.81 


-61: 


43:35.1 


18.91 


316.4 


21 c 


11:05:49.14 


-61: 


43:41.6 


19.59 


68.4 


22 


11:05:45.16 


-61: 


43:32.1 


19.76 


327.7 


23 


11:05:46.24 


-61: 


43:28.1 


20.93 


352.5 


24 


11:05:44.49 


-61: 


43:34.1 


21.31 


313.8 


25 


11:05:49.76 


-61: 


43:48.5 


22.65 


89.2 



a To convert to the 2MASS frame, add 0.17 s to a 
and add 071 to 5. 
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b With respect to WR 38 = target #1. 
C WR 38a. 
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Table 2. Uncorrected Photometric Data for the WR 38/WR 38a Cluster 



Number mF555W "T-F439W ~ m F555W ^F336W - "T-F439W V B - V 



l a 




15 


,00±0. 


,03 


1 


.21±0. 


.05 


A O T 1 A A O 

0.37±0.06 


14 


.97±0. 


03 


1 


,13±0.05 




O T^Tt A A 

blP40. 


















14.66 






1.28 


2 




16 


.47±0 


.05 


1 


.21±0. 


.13 


— 0.18±0.14 


16 


,44±0. 


06 


1 


,13±0.13 




STPH9 


















16.21 






1.23 


3 . . 




16 


.19±0 


.05 


1 


,09±0. 


1 A 

.10 


A AO I A 1 1 


16 


1 n i a 

.lo±0. 


a r 

.05 


1 


AO I A 1 A 

.03±0.10 




f~i run -i A-i 

bTPlOl 


















15.90 






1.27 


4 




17 


,60±0 


.12 


1 


.35±0 


.42 




17 


.57±0. 


.12 


1 


O A 1 A A O 

,24±0.42 


5 . . 




16 


,81±0. 


.06 


1 


.32±0 


.19 


0.13±0.24 


16 


.79±0, 


A^7 

.07 


1 


,21±0.19 




blP138 


















16.56 






1.33 


6 . . 




17 


.39±0 


.12 




















7 . . 
8 . . 




13 
18 


.27±0 
.79±0 


.03 

.39 


1 


,41±0 


.04 


A OA I A A A 

0.20±0.04 


13 


.25±0 


03 


1 


OA 1 A A r* 

,29±0.05 


9 . . 




18 


,70±0 


.40 




















10 




15 


.83±0 


.04 


1 


,21±0 


.09 


r\ r\c\ i r\ -i c\ 

0.09±0.10 


15 


.80±0, 


.04 


1 


.13±0.09 




blP76. 


















15.51 






1.28 


11 




18 


.63±0 


.32 




















12 




18 


.69±0 


.35 




















13 




16 


.96±0 


.07 





.85±0 


.16 


0.28±0.20 


16 


.94±0. 


.08 





,82±0.16 


14 




16 


.74±0 


.07 


1 


.42±0 


.22 


0.08±0.26 


16 


,72±0, 


07 


1 


.30±0.22 




STP137 


















16.51 






1.30 


15 




17 


,55±0 


.13 


1 


,41±0 


.43 




17 


,53±0. 


13 


1 


,29±0.43 


16 




17 


11±0 


.09 


1 


.24±0 


.25 


0.00±0.29 


17 


,09±0, 


09 


1 


,15±0.26 


17 




17 


.83±0 


.16 


1 


.22±0 


.47 




17 


,80±0. 


.16 


1 


,13±0.47 


18 




18 


,79±0 


.40 




















19 




16 


.00±0 


.04 


1 


.07±0 


.09 


-0.04±0.10 


15 


,98±0. 


.05 


1 


,01±0.09 


20 




17 


,99±0 


.17 





.83±0 


.39 




17 


,97±0. 


.17 





,80±0.39 


21 b 




15 


.33±0 


.04 


1 


.24±0 


.07 


-0.13±0.08 


15 


,30±0. 


.04 


1 


,15±0.07 




STP57. 


















15.12 






1.15 


22 




18 


12±0 


.21 




















23 




15 


.38±0 


.04 





,45±0 


.05 


0.42±0.05 


15 


,36±0. 


.04 





,45±0.06 


24 




18 


.22±0 


.23 




















25 




17 


.50±0 


.13 


1 


,23±0 


.38 




17 


,47±0. 


13 


1 


,15±0.38 



a WR 38. 



b WR 38a. 
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Note. — Corresponding stars and results from Shorlin ct al. (2004) (STP#) are listed below their 
HST values. 



- 21 - 



Table 3. ZAMS Intrinsic Colors from Synthetic Magnitudes 





m(336) - V 


m(439) - V 


m(555) - V 


U-V 


B-V 


J — V 


H-V 


K-V 


K 


10000 


-0.15 


-0.02 


-0.00 


-0.07 


-0.02 


0.06 


0.08 


0.08 


1.46 


15000 


-1.04 


-0.14 


-0.01 


-0.68 


-0.14 


0.34 


0.39 


0.44 


0.46 


20000 


-1.42 


-0.20 


-0.02 


-0.95 


-0.19 


0.49 


0.57 


0.63 


-0.44 


24000 


-1.63 


-0.22 


-0.02 


-1.10 


-0.22 


0.58 


0.67 


0.76 


-1.09 


27500 


-1.75 


-0.26 


-0.02 


-1.20 


-0.25 


0.63 


0.72 


0.82 


-1.66 


35000 


-1.90 


-0.28 


-0.02 


-1.32 


-0.27 


0.70 


0.81 


0.92 


-2.97 


45000 


-1.96 


-0.28 


-0.02 


-1.36 


-0.28 


0.71 


0.82 


0.93 


-4.79 
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Table 4. i? = 3.1 Reddening Parameters for Stars with 2MASS Photometry 



(m(336) - K) 



Field # 


2MASS # 


STP # 


E(B - V) 


-£(m(336) - K) 


K - A{K) 


Comment 


1 


11054642- 


6113487 


40 


1.46 ±0.12 


-1.35 ±0.06 


10.43 ±0.05 


WR 38 


3 


11054613- 


-6113457 


101 


1.51 ±0.09 


-2.68 ±0.10 


12.16 ±0.07 


all bands 


5 


11054738- 


6113550 


138 


1.70 ±0.09 


-2.81 ±0.17 


12.31 ±0.03 


all bands 


7 


11054719- 


6113398 




1.84 ±0.11 


-2.91 ±0.03 


8.36 ±0.02 




8 


11054805- 


6113469 




2.51 ±0.24 




11.97 ±0.06 




9 


11054664- 


6114003 




1.70 ±0.20 




14.24 ±0.12 




10 


11054806- 


6113526 


76 


1.61 ±0.10 


-2.71 ±0.09 


11.54 ±0.06 


all bands 


12 


11054562- 


6113358 




0.20 ±0.10 






foreground 


13 


11054515- 


6113367 




0.70 ±0.08 






foreground 


14 


11054845- 


6113414 


137 


1.69 ±0.15 


-2.62 ±0.17 


12.16 ±0.08 


all bands less J 


15 


11054869- 


6113465 




1.47 ±0.06 




13.21 ±0.09 




19 


11054659- 


6113301 




1.52 ±0.10 


-2.98 ±0.09 


12.13 ±0.06 




21 


11054896- 


6113413 


57 


1.49 ±0.09 


-1.60 ±0.06 


10.33 ±0.04 


WR 38a 


23 


11054607- 


6113279 




0.55 ±0.04 


-0.56 ±0.12 


13.84 ±0.11 


foreground 


25 


11054956- 


6113482 




1.39 ±0.09 




13.28 ±0.06 
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Table 5. Average Magnitude Residuals from Reddening Fits 



Filter < m A (observed) — m x (synthetic) > 



F336W 


0.20 ±0.05 


F439W 


0.04 ±0.04 


F555W 


0.23 ±0.04 


U 


-0.17 ±0.06 


B 


0.03 ±0.06 


V 


0.02 ±0.01 


J 


0.00 ±0.01 


H 


-0.06 ±0.09 


K 


-0.15 ±0.07 
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